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ABSTRACT 

This is the third paper in a series analyzing X-ray emission from the hot interstellar medium in a 
sample of 54 normal elliptical galaxies observed by Chandra. We focus on a subset of 36 galaxies with 
sufficient signal to compute radial temperature profiles. We distinguish four qualitatively different 
types of profile: positive gradient (outwardly rising), negative gradients (falling), quasi- isothermal 
(flat) and hybrid (falling at small radii and rising at larger radii). We measure the mean logarithmic 
temperature gradients in two radial regions: from 0-2 J-band effective radii Rj (excluding the central 
point source), and from 2-47? j. We find the outer gradient to be uncorrelated with intrinsic host 
galaxy properties, but strongly influenced by the environment: galaxies in low-density environments 
tend to show negative outer gradients, while those in high-density environments show positive outer 
gradients, suggesting the influence of circumgalactic hot gas. The inner temperature gradient, however, 
is unaffected by the environment but strongly correlated with intrinsic host galaxy characteristics: 
negative inner gradients are more common for smaller, optically faint, low radio-luminosity galaxies, 
whereas positive gradients are found in bright galaxies with stronger radio sources. There is no 
evidence for bimodality in the distribution of inner or outer gradients. We propose three scenarios to 
explain the inner temperature gradients: (1) Weak AGN heat the ISM locally, while higher- luminosity 
AGN heat the system globally through jets inflating cavities at larger radii; (2) The onset of negative 
inner gradients indicates a declining importance of AGN heating relative to other sources, such as 
compressional heating or supernovae; (3) The variety of temperature profiles are snapshots of different 
stages of a time-dependent flow, cyclically reversing the temperature gradient over time. 

Subject headings: galaxies: cooling flows — galaxies: elliptical and lenticular, cD — galaxies: ISM — X- 
rays: galaxies — X-rays: ISM 



1. INTRODUCTION 

In t he fir st two papers of this series (iDiehl fc Statlerl 
120071 120081 . hereafter Paper I and Paper II), we con- 
ducted a comprehensive morphological analysis of the 
hot gas in normal elliptical galaxies. In Paper I, we in- 
troduced a technique to separate the hot gas emission 
from the contamination of unresolved point sources. We 
applied this technique to a Chandra archive sample of 
54 elliptical galaxies and presented a gallery of adap- 
tively binned gas-only images, which were photon- flux 
calibrated and background corrected. We used these 
gas maps to derive isophotal cllipticity profiles and con- 
ducted a systematic morphological analysis. A compar- 
ison between optical and X-ray ellipticities measured in 
the inner, stellar mass dominated regions shows no cor- 
relation, contrary to what would be expected if the gas 
were in perfect hydrostatic equilibrium. We modeled 
the expected correlation under various assumptions, and 
concluded that these systems, in general, are at best only 
approximately hydrostatic. Moreover, the gas morpholo- 
gies almost always look disturbed. In Paper II, we in- 
troduced a quantitative measure of morphological asym- 
metry, and found it to be tightly correlated with radio 



Electronic address: diehl@lanl.gov, statler@ohio.edu 

1 Theoretical Astrophysics Group (T-6)/Computational Meth- 
ods Group (CCS-2), Mailstop B227, Los Alamos National Lab- 
oratory, P.O. Box 1663, Los Alamos, NM 87545, USA (present 
address) 

2 Astrophysical Institute, Department of Physics and Astron- 
omy, 251B Clippinger Research Laboratories, Ohio University, 
Athens, OH 45701, USA 



continuum power and nuclear X-ray luminosity. We also 
found the gas morphology to be influenced, to a compara- 
ble degree, by the ambient intergalactic medium. But the 
AGN-morphology correlation forms a continuous trend 
down to the lowest detectable AGN luminosities, indi- 
cating the importance of AGN feedback, even in rather 
X-ray faint elliptical galaxies. 

In this third paper, we address the question of 
whether the central AGN is merely redistributing 
the gas, or heating it as well. We produce ra- 
dial temperature profiles and find that they fall into 
a variety of distinct types. In particular, we con- 
firm th at negative (outwardly f a lling) temperatu r e gra- 
dients (iHumphrey et~aU 120061: iFukazawa et all 120061 : 
iRandall et al.ll2006UKhosroshahi et al.ll2004[ ) are present , 
and relatively common, in low- luminosity systems. 

Outwardly rising (positive-gradient) temperature pro- 
files, nearly ubiquitous in galaxy clusters, X-ray groups, 
and massive ellipticals, are usually understood as being 
the result of efficient radiative cooling in the dense cen- 
tral regions. Accreting gas at large radii can additionally 
shock- heat itself, amplifying the positive gradient. This 
interpretation is supported by the short central cooling 
times observed for galaxies, groups, and clusters, which 
can drop well below 100 Myr. However, cooling times are 
equally short in galaxies with ne gative temperature gra - 
dients, i.e. with a warm center (Humphr ey et al.1 f2006). 
Several solutions have been proposed to explain these 
c ounter-intuitive objects . 

Fukazawa et al. (2006) suggest that the gradients are a 
function of environment, with outwardly rising (positive) 
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gradients caused by the hotter ambient intracluster or 
intragroup gas surrounding the galaxy. Galaxies with 
n egative gradients s hould instead be isolated. 

iHumphrev et al.l (|2006( ), on the other hand, propose 
a bimodal distribution of temperature gradients, and 
suggest that the total mass of the system is the deci- 
sive factor for the sign of the gradient. The division 
between their two groups happens at a virial mass of 
~ 10 13 Mq, implying a distinction between normal galax- 
ies and groups. They hypothesize further that the tem- 
perature gradients could be related to a significant dif- 
fe rence in the galaxies' ey olutionary histories. 

Hum phrey et al.1 (2006) also discuss the role of com- 
pressive (gravitational) heating, noting that during a 
slow inflow of relatively cool gas (< 1 — 2keV) the en- 
ergy gain would exceed the radiative losses. For the in- 
flow of hotter baryons, radiative cooling would dominate 
and one would observe a positive temperature gradient 
instead. However, they find no reason for hot baryons to 
be specific to systems above their break mass ~ 10 13 Af Q , 
a nd suspect the envi r onme nt as a fuel source instead. 

iKhosroshahi et all (|2004f ) observe a negative tempera- 
ture gradient in the fossil group candidate NGC 6482 and 
argue along the same lines. They model NGC 6482's tem- 
perature profile successfully with a steady-state cooling 
flow with a reasonable cooling rate of M = 2M Q yr _1 , 
and adopt it as their preferred solution. They also es- 
timate that type la supernovae (SN) may be responsi- 
ble for balancing about 1/3 of the radiative losses in this 
galaxy. They find the contribution from type II SN to be 
insignificant and argue against AGN feedback on grounds 
of the very relaxed appearance of NGC 6482. 

In this paper, we show that the distribution of tem- 
perature gradients is not bimodal. We further show that 
the temperature gradients within the inner 2 optical ef- 
fective radii are not strongly influenced by the environ- 
ment. Instead, we find evidence that these inner gradi- 
ents owe their origins either to the specific nature of low- 
luminosity AGN feedback or to a declining importance 
of AGN relative to compressive heating or supernovae. 

In we summarize our analyses and results from Pa- 
pers I and II, and describe the methodology to derive 
radial temperature profiles. We then discuss the various 
types of temperature profiles seen in our sample in [J31 
For a quantitative analysis, we split the radial range into 
two regions: the inner region extending out to 2 effective 
radii and an outer region between 2 — 4 effective radii. 
We fit and analyze the gradients in these two regions 
separately, and demonstrate that the inner gradient is 
determined by galaxy properties, while the outer gradi- 
ent is strongly influenced by the presence of neighboring 
galaxies and/or a hot ambient medium. In |j4l we dis- 
cuss the implications of our findings for cooling flows, SN 
heating, and AGN feedback, before we briefly summarize 
in U3 

2. DATA ANALYSIS 

2.1. Summary of Results from Papers I and II 

We make use of several parameters from Papers I and 
II. We list those essential to our analysis in Tables [Hand 
[2] for completeness, along with some additional quanti- 
ties. We extract absolute K magnitudes Mk and J-band 
effective radii i?j from the 2MASS extended source cat- 



alog (jJarrett et al.l l2000h . We adopt 20 cm radio con- 
tinuum radio lumin osities -Lnvks from the NRAO VLA 
Sky Survey (NVSS. [Condon et al.|[l998h within 3i?, T (see 
Paper II) . In addition, we extract 6 cm radio continuum 
luminosities La cul fro m the GB6 catalog of radio sources 
(iGreeorv et all 11991. the Parkes-MIT-NRAO 4.85GHz 
Survey s (IWright et al.l I1996D . and a 6 cm radio cata- 
log by iBecker et al.l (|1991ft in the same region. Cen- 
tral velocity dispersion values are taken fr om the Lyon- 
Meud on Extragalactic Database (LEDA; iPaturel et al.l 
|1997|) . We also adopt the projected galaxy density pa- 
rameter P2MASS from Paper II, which is based on the 
number of neighbors in the 2MASS extended source cat- 
alog, and corrected for incompleteness. As it is one of 
the few accessible paramete rs to describ e the galaxy en- 
vironment, we also list the iTullyl (|1988h galaxy density 

PTully 

2.2. Temperature Profiles 

To produce radial temperature profiles, we divide the 
X-ray counts image of each galaxy into elliptical annuli, 
according to the X-ray ellipticity profiles computed in 
Paper I. For those galaxies with insufficient signal to fit 
ellipses, we revert to circular annuli. We find no evi- 
dence that this choice affects our results in any way. We 
adapt the width of our annuli to contain a minimum of 
900 counts above the background level, which we deter- 
mine by the appropriately rescaled Markevitch blank-sky 
background files 3 . 

We then extract a source and background spectrum for 
each annulus and fit them with a two-component model 
using the CIAO analysis package Sherpa. The first com- 
ponent consists of an APEC 4 plasma model to represent 
the hot gas emission. A quantitative comparison with 
its better-known predecessor, the Mckal model, shows 
nearly identical fitting results. We fix the gas metal- 
licity at the solar abundance value. Unresolved point 
sources are represented by a power-law model with the 
power law index fixed at 1.6. This "universal" spectral 
model is an adequate representation for the emission of 
low- luminosity low- mass X-ray binaries, as d emonstrated 
in Pap er I and determined independently bv llrwin et al.l 
( 2003). We also add a multiplicative absorption compo- 
nent, for which we fix the hydrogen column density to 
the Galactic value, evaluated at the target position with 
the CIAO tool C olden 5 . 

We repeat our spectral analysis for a few objects with 
the gas abundance as a free parameter, and find that 
our choice to fix them to the solar value does not affect 
the fitted temperature. Since the metallicity is poorly 
constrained by the fits in low signal-to-noise systems, we 
fix the metallicity for all of our galaxies, in order not to 
introduce systematic differences in the analysis. 

3. RESULTS 

3.1. Radial Temperature Profile Types 

We categorize the observed temperature profiles into 
four major groups, described below. Two examples from 
each are shown in Figure [TJ (Note that the distinctions 
between the groups are not always clear cut.) 

3 http:/ /cxc. harvard.edu/cal/Acis/Cal_prods/bkgrnd/acisbg/COOKBOOK 

4 Astrophysical Plasma Emission Code 

5 http:/ /cxc. harvard.edu/toolkit/colden.jsp 
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TABLE 1 

Chandra X-RAY LUMINOSITY AND TEMPERATURE PROFILE PARAMETERS. 
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TABLE 2 

Optical, radio, and environmental parameters 
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139 ± 2 


< 3.8 


X 


10 2b 


< 2.7 


X 


10 27 




0.49 


NGC3379 


-23.85 ± 0.11 


29.9 


205 ± 2 


3.2 ±0.7 


X 


10 26 


< 2.4 


X 


10 27 


3.82 ± 0.10 


0.52 


NGC3585 


-24.81 ± 0.18 


32.3 


207 ±4 


< 1.2 


X 


10 2/ 


< 2.0 


X 


10 28 


2.43 ± 0.31 


0.12 


NGC3923 


-25.30 ±0.28 


43.8 


247 ±6 


< 1.6 


X 


10 27 


< 2.6 


X 


10 28 


2.89 ± 0.16 


0.40 


NGC4125 


-25.03 ± 0.25 


33.0 


226 ±6 


1.7 ±0.4 


X 


10 28 


< 1.2 


X 


10 28 


3.13 ± 0.13 


0.34 


NGC4261 


-25.24 ±0.19 


25.5 


320 ±8 


1.0 ±0.2 


X 


10 31 


4.8 ±0.1 


X 


10 30 


3.09 ± 0.14 


0.84 


NGC4365 


-24.91 ± 0.17 


40.7 


255 ± 2 


< 1.2 


X 


10 27 


< 9.0 


X 


10 27 


3.16 ± 0.13 


2.93 


NGC4374 


-25.10 ±0.11 


34.8 


280 ± 2 


2.5 ±0.3 


X 




1.3 ±0.0 


X 


10 30 


3.18 ± 0.14 


3.99 


NGC4406 


-25.07 ±0.14 


59.7 


235 ± 2 


1.1 ±0.2 


X 


10 27 


< 6.3 


X 


10 27 


3.31 ± 0.13 


1.41 


NGC4472 


-25.66 ± 0.10 


59.2 


288 ± 2 


8.1 ±0.8 


X 


10 28 


2.7 ±0.0 


X 


10 28 


3.41 ± 0.12 


3.31 


NGC4494 


-24.16 ±0.11 


30.8 


149 ±3 


< 8.7 


X 




< 6.3 


X 


10 27 




1.04 


NGC4526 


-24.67 ±0.20 


43.8 


263 ± 18 


6.6 ± 1.2 


X 


10 27 


< 6.1 


X 


10 27 


2.54 ± 0.31 


2.45 


NGC4552 


-24.20 ±0.14 


25.4 


253 ±2 


2.9 ±0.4 


X 


10 2S 


1.9 ±0.0 


X 


10 28 


3.62 ± 0.09 


2.97 


NGC4555 


-25.78 ± 0.33 


10.9 




< 2.8 


X 


10 28 


< 2.0 


X 


10 29 


3.19 ± 0.22 




NGC4564 


-22.94 ±0.17 


19.9 


158 ± 2 


< 6.7 


X 


10 26 


< 4.8 


X 


10 27 




4.09 


NGC4621 


-24.56 ± 0.20 


32.9 


225 ±3 


< 10.0 


X 


10 26 


< 7.2 


X 


10 27 


2.66 ±0.25 


2.60 


NGC4636 


-24.41 ± 0.13 


59.3 


202 ±3 


2.7 ±0.3 


X 


10 28 


1.8 ±0.0 


X 


10 28 


3.47 ±0.12 


1.33 


NGC4649 


-25.39 ± 0.15 


45.2 


334 ±3 


9.8 ± 1.4 


X 


10 27 


1.3 ±0.0 


X 


10 28 


3.39 ±0.12 


3.49 


NGC4697 


-23.98 ± 0.14 


42.4 


173 ± 2 


< 4.1 


X 


10 26 


< 6.6 


X 


10 27 


3.65 ±0.11 


0.60 


NGC5018 


-25.27 ±0.33 


15.6 


214 ±8 


< 4.8 


X 


10 27 


< 8.0 


X 


10 28 


2.50 ±0.31 


0.29 


NGC5044 


-24.76 ± 0.28 


25.3 


238 ±8 


4.0 ± 1.0 


X 


10 28 


< 4.9 


X 


10 28 


3.17 ±0.13 


0.38 


NGC5102 


-21.09 ± 0.14 


79.3 


66 ±4 


6.1 ± 1.5 


X 


10 25 


< 1.4 


X 


10 27 




0.17 


NGC5171 


-24.95 ± 0.33 


10.8 




< 2.9 


X 


10 28 


< 2.1 


X 


10 29 






NGC5532 


-26.33 ± 0.33 


16.0 


293 ± 18 


5.8 ± 1.7 


X 


10 31 


1.6 ±0.1 


X 


10 31 


3.11 ±0.25 




NGC5845 


-22.96 ± 0.21 


4.9 


234 ±8 


< 2.0 


X 


10 27 


< 1.4 


X 


10 28 




0.84 


NGC5846 


-25.04 ±0.20 


34.5 


239 ±3 


1.6 ±0.3 


X 


10 28 


< 1.3 


X 


10 28 


2.97 ±0.15 


0.84 


NGC6482 


-25.48 ± 0.33 


12.6 


303 ±9 


< 1.0 


X 


10 28 


< 7.5 


X 


10 28 


2.64 ±0.31 




NGC7052 


-25.66 ± 0.33 


21.8 


271 ±9 


1.2 ±0.4 


X 


10 30 


6.7 ±0.6 


X 


10 29 


2.72 ±0.31 




NGC7618 


-25.40 ±0.33 


11.6 




2.7 ±0.8 


X 


10 29 


< 1.3 


X 


10 29 


2.76 ±0.31 





a 2MASS data from the extended source catalog ( Jarrctt ct al. 2000). _R"-band absolute magnitude Mk and J-band 
effective radius i?j in kpc. b Velocity dispersion (in km s" 1 ) from LEDA (Paturcl ct al. 1997).° LjmvssO 11 er S s 
20cm radio continuum luminosity from NVSS (Condon ct al. 1998); Lecm (in erg s _1 ): Combination of 6cm radio 
continuum luminosities from Gregory et al. ( 1996), Wright ct al. ( 1996), and Becker et al. ( 1991). d Projected galaxy 
density P2MASS i s scaled logarithmically and in units of Mpc — 2 , Tully galaxy density pxully i s scaled linearly and in 
units of Mpc - 3 . 
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Positive Gradients. — These temperature profiles show a 
positive gradient at all radii, i.e. the temperature con- 
tinuously rises outwardly. These profiles resemble those 
found in clusters of galaxies, which generally harbor cool 
cores. 

Negative Gradients. — The temperature profiles show 
a negative gradient at all radii, i.e. tempera- 
tures monotonically decline outward. This phe- 
nomenon is l ess well-known, and has been reported 
only recently dRandall et al.l 2006; Fukaz awa et alj|2006t 
Humphrey et alj|2006t iKhosroshahi et al.ll2004h . 

Hybrid. — These peculiar cases exhibit a dramatic 
change in the temperature gradients. The gradient 
changes its sign from negative to positive at some inter- 
mediate radius, generally between 1 — 3i?j. These galax- 
ies have warm centers, outside of which their tempera- 
tures drop to a minimum level and ris e back up again. 
These profiles have first been noted bv lHumphrev et"aT1 
((200(1 . 

Quasi-Isothermal. — The radial temperature profiles in 
this category are consistent with being almost flat at all 
radii. These galaxies form the transition point between 
galaxies with positive and negative gradients. 

We observe only hybrid temperature profiles that 
change their gradient from negative to positive. Some 
cooling flow clus ters have been found to exhibit the op- 
posite behavior (|Piffaretti et alj |2005|) . Their temper- 
ature profiles show a cool center, then rise to a peak 
temperature and fall back down on the outskirts. This 
"break" usually happens at around 10% of the virial ra- 
dius, which is larger than the radii tha t we are probing 
in nor mal galaxies. A ROSAT study bv I O 'Sullivan et all 
(2003) exhibits similar trends for elliptical galaxies at 
larger radii. 

We split the profiles into two radial regions and analyze 
the inner and outer temperature gradients separately. As 
most hybrid profiles exhibit their turnover in slope some- 
where around 2i?j, we use this radius as the boundary 
between our two regions. Accordingly, we define the in- 
ner region from outside the central point source extend- 
ing out to 2i?j and the outer region between 2 — 4i?j. 
We then fit each part of the profile with a power law 
to derive effective temperature gradients for each region. 
We will refer to the logarithmic gradients dlnT/dlni? 
evaluated within 2Rj and from 2 — 4i?j as «02 and «24, 
respectively. 

The best-fit values for 0:02 and «24 are listed in Table 
[TJ The reported errors are the formal la statistical er- 
rors obtained from the fitting procedure. For cases with 
only 2 valid data points within the fitting range, we use 
the difference between these two points to derive a gradi- 
ent; the errors are derived by propagating the statistical 
errors of the individual temperature measurements. 

3.2. The Inner Temperature Gradient «02 

Figure [2] shows the compilation of all profiles overlaid 
in one plot, with the radial axes scaled by their J-band 
effective radius Rj. This plot already clearly indicates 
the absence of any real bimodality in the temperature 
profiles. Each profile is colored according to its lumi- 
nosity weighted temperature Tx within 3 optical radii. 
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Fig. 1. — Examples of different projected temperature profiles 
as a function of radius. Temperature profile types can be divided 
into 4 major groups (top to bottom rows): (1) Positive gradient 
(outwardly rising) at all radii; (2) Negative gradient (outwardly 
falling) at all radii; (3) Hybrid, negative gradient in the core and 
positive gradient at larger radii; (4) Quasi-isothermal, no apparent 
temperature change with radius. The complete set of temperature 
profiles is available online for all 36 galaxies with two or more valid 
temperature profile points. 

The coloring changes smoothly from top to bottom, in- 
dicating overall positive gradients for intrinsically hotter 
galaxies and negative gradients for cooler galaxies. We 
plot ao2 as a function of Tx in Figure [3J Again, there 
is no sign of bimodality in ao2- The plot shows a tight 
correlation, significant at the 9.1a level, reflecting the 
fact that we observe only a very small range in central 
temperatures between 0.6 — 0.7keV (Figure [2]). Thus, 
any average temperature will obviously be strongly cor- 
related with the gradient as well. Consistently, the fit 
suggests that the transition from positive to negative in- 
ner temperatures gradient occurs at a mean temperature 
around 0.64 keV. 

To establish the underlying cause for the negative inner 
temperature gradients, we perform a correlation analy- 
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Fig. 2. — Combined plot of all projected temperature profiles, 
as a function of radius scaled by the J- band effective radius Rj. 
Error bars are omitted for clarity; for typical error estimates, refer 
to Figure [T] The profiles are colored according to the luminosity 
weighted average temperature of the galaxy within 3 optical radii 
Tx, as indicated by the color scale bar. The temperature gradient 
changes continuously from positive gradients at the top to isother- 
mal and hybrid profiles to negative gradient profiles at the bottom, 
along with the average temperature. 

sis with various galaxy properties, the most interesting 
of which are listed in the upper half of Table [3J X-ray 
gas luminosities £x,gas, absolute K magnitudes Mk, cen- 
tral velocity dispersions a, radio luminosities at 20cm 
(Lnvss) and 6cm (L§ cm ), and environmental measures 
of local galaxy density p 2 MASS and pnaiy. 

We assess the correlation of «02 with each of these 
properties using the linear fitting algorithm bandf it (see 
Appendix of Paper II). This algorithm models the dis- 
tribution of (x, y) points as a linear band with a Gaus- 
sian intrinsic width. The model is fitted by maximiz- 
ing the likelihood of the data; this is similar to fitting a 
straight line by minimizing the error-weighted perpendic- 
ular residuals. The standard error a in the best-fit slope 
is obtained from the covariance matrix, and the signif- 
icance of the correlation is the number of a by which 
the slope differs from zero. In these fits, the quantities 
assigned as abscissae are scaled logarithmically (except 
for Mk, which is already intrinsically logarithmic), while 
the ordinate an2 is scaled linearly. The top half of Ta- 
ble [3] lists the best-fit parameters and the statistical sig- 
nificance of the correlation. The fitted parameters are 
omitted where the significance is < 2a. The tabulated 
fits are obtained from the full data set including upper 
limits, but we also include in the last column the signif- 
icance of the correlation obtained with the upper limits 
omitted. 6 The xq column indicates the transition point 
from negative to positive temperature gradients in each 
fit. 

As the table and Figure 0] show, the strongest corre- 
lations with ao2 are found with the 20 cm NVSS radio 
luminosity £nvsSj the- velocity dispersion er, and the ab- 
solute K-magnitude Mk- All of these three correlations 
are of roughly the same significance, with comparable in- 
trinsic widths, as Figure [4] shows. The correlation with 

6 The exception is Lgb6; f° r which we give the parameters for 
the fit without upper limits. This is because the number of galaxies 
with upper limits only (49) greatly outweighs that of galaxies with 
actual detections (15). 



the 6 cm radio luminosity Lgb6 seems equally tight, but 
has a lower significance, due to the smaller sample size 
with 6 cm radio luminosity measurements. 

We do not find any evidence that ao2 is correlated 
with the environmental galaxy densities /92MASS or pTully 
The /32MASS _ «02 plot is shown in the bottom panel of 
Figure O We conclude that the inner temperature gra- 
dients are neither the result of interactions with neigh- 
bor galaxies, nor with ambient intragroup or intraclus- 
ter gas. Instead, we find that they are connected to 
intrinsic galaxy properties. We can generally character- 
ize galaxies with negative inner temperature gradients as 
being smaller, optically fainter galaxies with lower veloc- 
ity dispersions, lower X-ray gas luminosities, lower aver- 
age temperatures, and lower radio luminosities than their 
positive gradient counterparts. 

Unfortunately, all of these galaxy properties are 
intimately connected with each other through well- 
known correlati o ns su ch as the Tx~a relation (e.g. 
IP 'Sullivan et all f2003h the Faber- Jackson relation 
(Faber & Jac ksqnl Il976f) . the £x, gas -Tx relation (e.g. 
IQ 'Sullivan et al. [|2g03D and the I'Radio - c relation (e.g. 
iSnellen et alT feoO.U This makes it difficult to distin- 
guish between fundamental correlations that are really 
responsible for determining the inner temperature struc- 
ture and others that are simply "riding along" via other 
correlations. 

We check the robustness of our results by deriving in- 
ner temperature gradients for different cutoff radii and 
find that all of the observed trends are confirmed, as long 
as the cutoff-radius does not exceed ~ 3i?j. In particu- 
lar, we find that for smaller cutoff-radii (e.g. li?j), the 
significance of the correlations with a, Mk and £x,gas 
slightly decreases, while the correlations with the radio 
luminosities Lnvss and Lq cu1 strengthen even further. 
This may suggest that the correlations with radio lumi- 
nosities are intrinsically the strongest. Figure [6] shows a 
combined plot of all temperature profiles, similar to Fig- 
ure but this time colored according to the NVSS radio 
luminosities. A trend with radio luminosity is clearly ev- 
ident. We will discuss the implications of our results in 

SI 

3.3. The Outer Temperature Gradient «24 

We now look at the outer temperature gradient be- 
tween 2 and 4 Rj , and with what it is correlated. Like 
<*02j Q!24 is correlated with the average temperature 
within 3-Rj (Fig. [7]), though less strongly (2.5a, com- 
pared with 9.1a). This trend is such that galaxies with a 
hotter average temperature have stronger positive tem- 
perature gradients. Some correlation is expected because 
these quantities are not truly independent. As a conse- 
quence, the outer temperature gradient is also weakly 
correlated (2.4ct) with the inner temperature gradient. 

We now repeat the same analysis for the outer temper- 
ature gradient 024 . The results of the correlation analysis 
are listed in the bottom half of Table [3] Unlike ao2 , «24 
does not depend on the intrinsic galaxy properties Lnvss, 
Lecm, c, or Mk- Instead, we find strong evidence that 
«24 depends only on the environmental density param- 
eters P2MASS and pTuiiy This trend with environment 
is statistically even stronger than the one with the lu- 
minosity weighted temperature Tx, even though those 
parameters are not independent measurements. 
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TABLE 3 

Correlations involving inner and outer temperature gradients. 



y 




(2 





XQ 


Significance 
(with limits) 


Significance 
(without limits) 




log Lwv^ 1 ; 


0.060 ± 0.019 


-1.65 ± 0.02 


27 .4 


3.2ct 


4.3ct 




log <T 


4.171 ± 1.192 


-10.0 ±0.01 


2.40 


3.6<t 


3.6cr 


Q:()2 


M K 


-0.141 ± 0.044 


-3.50 ±0.02 


-24.8 


3.2a 


3.2cr 


«02 


logL G B6* 


0.083 ± 0.036 


-2.35 ±0.04 


28.3 


0.4ct 


2.3cr 




lo£T /jY 








I.Oct 


I.Oct 


«02 


logp 








0.7ct 


0.7fJ 


«02 


log PTully 








0.4ct 


0.4er 


Q24 


logp 


0.735 ± 0.223 


-2.28 ±0.03 


3.10 


4.4ct 


4.4ct 


024 


log PTully 


0.272 ± 0.102 


0.04 ± 0.04 


0.15 


2.8ct 


2.8er 


«24 


M K 








1.6CT 


1.6(7 


Q24 


logo- 








0.5ct 


0.5cr 


C*24 


log L GB6 








< O.lcr 


0.6cr 


024 


log I/NVSS 








< O.lcr 


0.3fJ 


024 


log £x,gas 








< O.lcr 


< O.lcr 



Note. — Results are listed in order of decreasing correlation significance. Parameters refer 
to linear fits of the form y = ax + b (for correlations of > 2cr significance), xo denotes the 
point where the fit yields 0, i.e. where the temperature gradients change sign. The last two 
columns quote the correlation significances with and without upper and lower limits included, 
respectively. 

The fitted parameters for the -Lgb6 — a 02 correlation are quoted for the fit without upper 
and lower limits included. 

We conclude that the inner and outer temperature gra- 
dients are essentially decoupled. While the inner gradi- 
ent depends only on intrinsic galaxy properties, the outer 
gradient shows no correlations but with the environment. 

4. DISCUSSION 

4.1. Implications for Cooling Flows 

Steady-state cooling flow models have gone out of fash- 
ion recently due to extensive work on galaxy clusters, 
which show insufficient amount s of cooling gas at the 
center (|Peterson fc Fabian] [2006) . These simple models 
are unlikely to apply to X-ray bright elliptical galaxies 
either. We showed in Paper II that the hot gas in these 
systems is almost always disturbed, and we see evidence 
linking the origin of these disturbances to the central 
AGN. 

However, it is far from proven that the same is true for 
low-luminosity galaxies, in which we find negative inner 
temperature gradients. Compressive heating during a 
gradually cooling inflow of relatively cool gas may be 
able to offset radiative losses for l ow-temperature gas in 
steep gravitational potentials (e.g. lMathews fc Brighentl 
2003). This counter- intuitively results in a cooling flow 
that gets heated during inflow and may even produce 
a hot center, i.e. a negative gradient. We also find in 
Paper II that these systems are generally less disturbed, 
which could be consistent with a steady state cooling flow 

s olution. 

iKhosroshahi et ail (|2004D observe a falling tempera- 
ture profile for the fossil group candidate NGC 6482 and 
successfully fit a steady-state cooling flow model with a 
reasonable cooling rate of M = 2M Q yr _1 . However, 
they derive the inner gravitational potential from the 
X-ray profiles themselves assuming hydrostatic equilib- 
rium, which yields a steep inner potential gradient. They 
then use this potential to fit the cooling flow model to 
the negative temperature gradient. This could be circu- 
lar reasoning: a central temperature peak together with 
the assumption of hydrostatic equilibrium implies a steep 
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Fig. 3. — Inner temperature gradient within 2 Rj as a function 
of the average luminosity weighted temperature within 3 optical 
radii. The dashed line indicates the best fit. 

To check the robustness of these results, we repeat our 
analysis using larger outer radial boundaries and confirm 
all trends. The significance of the environmental depen- 
dence gets even stronger when extending the analysis to 
larger radii. These relations are strongest, when one fits 
temperature gradients to all radii beyond 2i?jwithout 
imposing an outer radial limit. However, since the gra- 
dients tend to get stronger with radius, and our galaxies 
have very different cutoff radii owing to different surface 
brightness profiles, we do not report the functional form 
of the fit, as it is driven by the brightest galaxies. Nev- 
ertheless, this strengthens the confidence in the observed 
correlation. 
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Fig. 4. — Inner temperature gradient within 2 R j as a function of 20 cm NVSS radio luminosity (top left), central velocity dispersion 
(top right), absolute K magnitude (or stellar mass, bottom left) and 6cm radio luminosity (bottom right). The dashed lines indicate the 
best-fit correlations from bandfit, as given in Table [3] 



gravitational potential, which leads to increased com- 
pressional heating (e.g. [Mathews fe Brighentill2003l ) and 
a central temperature peak. In addition, we have argued 
in Paper I that an assumption of hydrostatic equilib- 
rium is generally not secure. It would be much safer to 
derive the gravitational potential from independent stel- 
lar dynamics, as the inner region is most likely stellar- 
mass dominated, at le ast within two effective radii (e.g. 
iMamon fe Lokasjl2005l ). In any event, NGC 6482 is the 
only negative-gradient object that has been successfully 
fitted with a cooling flow model so far. Only modeling a 
more complete sample will show if this idea can generally 
hold. 

4.2. Implications for the Existence of Circumgalactic 

Gas 

Negative tem perature gradients have been recognized 
only recently ( Humphrey et al. 2006; Fukazaw a et al.l 
l2006HRandall et alJl2006t iKhosroshahi et al.ll200l . Be- 
cause earlier observations revealed only positive gradi- 
ents, many theoretical flow models have been dismissed 



on the grounds that they pro duce negative gradients (e.g. 
iMathews fc Br ighenti 2003). Instead, theoretical effort 
has focused on findin g an explanation for the prev alence 
of positive gradients. [Brighenti fc Mathew s (1998) argue 
that a hot circumgalactic gas reservoir is able to reverse a 
negative temperature gradient. This explanation is con- 
sistent with our observations that the outer temperature 
gradient is correlated with the environment. Whether 
models with circumgalactic gas can quantitatively ac- 
count for the more complex hybrid temperature profiles 
remains to be seen. 

4.3. Implications for Supernova Feedback 

Another means of producing negative temperature gra- 
dients involves supernova (SN) feedback. Since star for- 
mation should be negligible in elliptical galaxies, this 
mechanism would involve only contributions from type 
la SN . Early proposed wind models involving SN feed- 
back (jBinnev fc Taborlll995() were later dismissed, since 
a main feature was a negative temperature gradient 
throughout their evolution, which had not been observed 
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Fig. 5. — Projected galaxy number density P2MASS vs - outer 
(c«24, top panel) and inner (oo2i bottom panel) temperature gra- 
dients. Note that ao2 is evidently unaffected by environment, 
whereas o?24 depends strongly on the density of nearby systems, 
suggesting the influence of hot ambient gas. 
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Fig. 6. — Temperature profiles, as in Figure [2] but with colors 
indicating NVSS 20 cm continuum radio luminosity within 3i?.j, as 
shown by the color bar. The radio luminosity changes continuously 
as profiles change from positive to negative temperature gradients. 
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Fig. 7. — Outer temperature gradient within 2 — 4i?j as a func- 
tion of the average luminosity weighted temperature within 3 -Rj . 
The dashed line indicates the best fit. 



at that time. However, iMathews fc Brighentil (|2003l ) 
point out that these models are able to reproduce ob- 
served gas profiles only for a very short period of time 
(~ 10 s yr) just before a cooling catastrophe sets in. Fur- 
thermore, these models are sensitive to the assumed SN 
rate, resulting in abrupt transitions to SN driven winds, 
and thus require fine-tuning. This fine-tuning problem 
can be circumvented by the pre sence of circumgalactic 
gas (|Mathews fc Brighen ti 2003). However, our observa- 
tions of purely negative gradients without an outer rise 
in temperature may present a challenge to this model, 
though we cannot exclude the possibility of an external 
gas reservoir below our detection limit. 

In addition, the predicted metallicities in galactic wind 
models are generally super-solar, significantly exceeding 
the historicall y observed extreme ly low abundances in 
the hot gas (|Arimoto et al.lll997f ). A Chandra spectral 
analysis of abu ndance gradients in a samp le of 28 ellipti- 
cal galaxies by lHumphrev fc Buotd (|2006l ). on the other 
hand, no longer shows strongly sub-solar abundances. 
They attribute this difference to previously imperfect 
modeling of the spectra, mainly caused by the neglect of 
the unresolved point source component and attempting 
to fit multi-temperature gas with a single-tem perature 
model, the so-called iron-bias (e.g. lBuotgl2000l ). Never- 
theless, they conclude that their abundances are still far 
too low to be consistent with gal actic wind mode l s, and 
favor the circulation flow model of lMathews et all (2004) 
instead. 

However, if we consider only the energy input 
by SNIa feedback, we find that SN feedback could 
play a role in heating the hot gas. The average 
SNIa rate for elliptical and SO gal axies is trn = 
0.18 ± O.Oe(lOOyr)- 1 (IO^Lbq)- 1 dCappellaro et all 
fl99flh . With an average energy injection of ~ 10 ergs 
per SN, this results in a SN heating rate of Lsn = 5.7 x 
10 30 (Lb/£b,©) ergs -1 . An inspection of the Lx, gas -LB 
diagram in Paper I shows that Lsn > Lx, gas for galaxies 
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below the blue luminosity where the inner temperature 
gradients turn from negative to positive (~ 4x 10 10 Lb©)- 
Although this is somewhat suggestive, the large scatter 
in X-ray luminosity of almost 2 orders of magnitude at 
a given blue luminosity, combined with the also rather 
large scatter in the Lb-Q.q2 relation, make it impossible 
to tell if this is simply coincidence. 

However, if SN heating were the main cause, we 
would expect a correlation between the inner temper- 
ature gradient and the SN heating to X-ray cooling ratio 
(XsN/£x,gas), in the sense that negative gradients would 
correspond to high ratios of heating to cooling. We do 
not find such a correlation. £sN/-^x,gas correlates with 
the inner temperature gradient very weakly on the less 
than 0.5a level. We conclude that supernova feedback 
may be important for balancing part of the radiative 
losses in X-ray faint galaxies, but our analysis suggests 
that it is most likely not the dominant factor. 

4.4. Implications for AGN feedback 

In Paper II we have measured the amount of asymme- 
try in the hot gas, and find a strong correlation between 
asymmetry and AGN power. This correlation persists 
all the way down to the weakest AGN luminosities at 
the detection limit of the NVSS 20 cm survey. We now 
find that the temperature structure is also strongly cor- 
related with the AGN luminosities, another hint toward 
the importance of AGN feedback throughout the ellipti- 
cal galaxy population. However, we cannot completely 
rule out compressionally heated cooling flows or SN feed- 
back to explain the prevalence of negative temperature 
gradients. Thus, we propose three possible scenarios in- 
volving AGN feedback to explain our results: 

1. Weak AGN with smaller black holes heat the ISM 
locally, while higher- luminosity sources feed power- 
ful jets that distribute the heat globally by blowing 
large cavities into the ISM. This is consistent with 
the observation that smaller elliptical galaxies have 
rather weak AGN and generally less extended radio 
emission, and also in agreement with our findings 
from Papers I and II that the amount of asymme- 
try correlates with AGN luminosity. In this sce- 
nario, weak AGN would still be disturbing the gas, 
but on a scale and surface brightness level that is 
simply less detectable, resulting in a lower asym- 
metry. Negative temperature gradients could then 
be a sign of very localized heating by the central 
AGN. 

2. AGN are responsible for globally heating the hot 
gas only in X-ray bright galaxies with positive tem- 
perature gradients. The onset of negative inner 
temperature gradients marks the point where AGN 
heating becomes unimportant, relative to other 
sources. These other sources could include com- 
pressional heating or supernovae. 

3. The observed temperature gradients are snapshots 
of different stages of a time-dependent flow, which 
cyclically reverses the temperature gradient over 
time. If such solutions exist, it will be challeng- 
ing for theoretical models to explain the fact that 
none of our galaxies exhibit central temperatures 



below ~ 0.6 keV. Thus, any cyclic solution has 
to keep the central temperature rather constant, 
while reversing the temperature gradient by heat- 
ing or cooling only at large radii. The best chance 
to achieve this may be for galaxies to cycle through 
wind and inflow phases, possibly intimately corre- 
lated with the time-dependent AGN activity of the 
central black hole. 

The possible importance of AGN heating for elliptical 
galaxi es has also recently been pointed out bv lBest et al.l 
(2006). They combine two empirical results to derive 
an estimate of time-averaged he ating by radio sourc es in 
galaxies. They use a result by iBirzan et all (|2004t l for 
galaxy clusters that empirically links the pdV work asso- 
ciated with inflating X-ray cavities into the intracluster 
medium with the observed 20 cm radio continuum power 
of the associated radio source. Although this correlation 
exhibits significant scatter, Best et al. derive a linear fit 
and use it to convert their radio powe r s to m echanical 
energy. In an earlier study, iBest et ail (|2005f ) find that 
the fraction of elliptical galaxies hosting radio-loud AGN 
correlates with black hole mass and radio luminosity. As- 
suming that all elliptical galaxies have AGN at their cen- 
ters, Best et al. interpret the fraction of galaxies with 
active AGN as the fraction of time that they are turned 
on. By combining the computed mechanical work per 
unit radio luminosity derived from IBirzan et al.l (2004) 
with the fraction of time the radio source is turned on, 
Best et al. calculate the time-averaged mechanical heat 
input of the AGN. A comparison with the Lx, s & s -Lb rela- 
tion for normal ellipticals shows a remarkable agreement 
between the time-average AGN heat input a nd the av- 
erage d radiative losses of elliptical galaxies (Best et al. 
2006, Figure 2). This good agreement is actually sur- 
prising, since the conversion factor from radio power to 
mechanical energy has a rather large scatter and is only 
based on observations of cluster cavities. 

Further support fo r AGN heating has been provided 
bv lAllen et all (|2006h . who measure the mechanical en- 
ergy associated with X-ray cavities in 9 X-ray luminous 
elliptical galaxies. They compare this value to the Bondi 
accretion rate, which they derive from deprojected den- 
sity and temperature profiles, evaluated at the accretion 
radius. Allen et al. find a tight correlation between the 
Bondi accretion rate and the mechanical energy injected 
into the ISM, and find that this energy input may be 
sufficient to prevent the gas from cooling. 

4.5. What is so special about ~ 0.6 keV? 

A close inspection of Figure[2]shows a remarkably small 
range in central temperature, which falls between 0.6 to 
0.7 keV. The upper limit owes its origin to our explicit 
exclusion of brightest cluster galaxies, with higher tem- 
peratures, from our sample. Including cluster cDs in our 
sample would add the missing profiles, adding positive 
temperature gradients with higher central temperatures. 

However, the lower limit is quite mysterious. We find 
it unlikely that this is simply a Chandra sensitivity ef- 
fect. We know that our temperature fits are sensitive 
to lower temperatures, as we can see them in fits to 
the outer regions of the same objects. Conceivably, this 
could represent a selection effect imposed on the Chan- 
dra archive through the proposal process, which disfavors 
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observations of systems with lower temperature due to 
the drop in instrument sensitivity at lower energy. We 
find this explanation also difficult to believe, as galax- 
ies with negative gradients would have been character- 
ized simply as having a lower mean temperature, since 
ROSAT would not have been able to detect the rise in 
temperature toward the center. However, we do see that 
the faintest galaxies in our sample exclusively build the 
lower envelope in the temperature profiles, with luminos- 
ity weighted temperatures of ~ 0.4 keV. Thus, fainter 
galaxies could lower this envelope even further, and with 
it the central temperature. The lower envelope may also 
mark the transition to a galactic wind, which would ren- 
der the temperature gradient for these galaxies unobserv- 
able due to low gas densities. 

Nevertheless, something is special about ~ 0.6 keV. 
First, we do not observe any central temperature below 
this value. Second, all hybrid temperature profiles drop 
below 0.6 keV at some intermediate radius and then rise 
back up again. And third, the best fit for the Tx-ao2 
relation puts the transition between negative and positive 
gradients at 0.64 keV. Any flow model on the galaxy scale 
has to be able to reproduce these properties. 

5. CONCLUSIONS 

We have reported on the shape of temperature pro- 
files in 36 normal elliptical galaxies. These profiles show 
a variety of different profile types: purely positive gra- 
dients, purely negative gradients, quasi- isothermal and 
even hybrid profiles. To understand this complexity, we 
derive mean temperature gradients for an inner region 
within 2i?j, excluding the central point source, and an 
outer region between 2 — 4i?j. We find that the outer 
temperature gradient is independent of intrinsic galaxy 
properties, but a strong function of environment, such 
that positive outer temperature gradients are restricted 
to cluster and group environments. This suggests that 
the outer gradients are caused by interaction with hot- 
ter ambient gas, whereas galaxies with negative outer 
gradients are in less dense environments and lack this 
intergalactic gas reservoir. 

The inner temperature gradient, on the other hand, is 
completely independent of the environmental influence. 
Instead, we find that it is correlated with a number of 
intrinsic galaxy properties; in decreasing order of signif- 
icance, the 20 cm radio luminosity, the central velocity 
dispersion, the absolute K magnitude, and the 6 cm ra- 
dio luminosity. 
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The data cannot rule out the idea that negative gra- 
dients can be produced by compressional heating in low- 
temperature systems, during a slow cooling inflow in 
a steep gravitational potential. SN feedback may also 
provide sufficient energy to offset cooling in X-ray faint 
galaxies, but we find no direct evidence that SN heating 
dominates. 

Our preferred feedback model involves the central 
AGN. The inner temperature gradient is most strongly 
correlated with radio luminositiy and central velocity dis- 
persion, which may be interpreted as a surrogate for 
black hole mass (|Tremaine et al.l l2002f ) . The nature of 
these correlations is such that weak AGN hosts show 
negative temperature gradients, whereas more luminous 
AGN exclusively live in positive gradient systems. Thus, 
we propose three scenarios, to explain the observed fea- 
tures. (1) Weak AGN distribute their heat locally, 
whereas luminous AGN heat the gas more globally with 
their extended jets. (2) The onset of negative gradients 
marks the point where AGN heating becomes unimpor- 
tant, and compressional heating or SN feedback becomes 
dominant. (3) A cyclic model in which the AGN drives 
an outflow, which shuts the AGN activity off until the 
flow reverses itself, fuels the black hole and starts another 
cycle. 

These findings are in agreement with the results from 
Paper I, which showed that precise hydrostatic equilib- 
rium does not hold for the hot gas in elliptical galax- 
ies, and established the prevalence of disturbances in 
the X-ray gas morphology. The results of Paper II in- 
dicate that the central AGN probably causes these dis- 
turbances. Combining these results with the connection 
between the temperature structure and the radio lumi- 
nosity of the system produces a strong argument for the 
general importance of AGN feedback in nearly all normal 
elliptical galaxies. 
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APPENDIX 
INDIVIDUAL TEMPERATURE PROFILES 
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